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SCHRODINGER EQUATIONS WITH PERIODIC POTENTIALS 
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Abstract. We consider semiclassically scaled Schrodinger equations with an 
external potential and a highly oscillatory periodic potential. We construct as- 
ymptotic solutions in the form of semiclassical wave packets. These solutions 
are concentrated (both, in space and in frequency) around the effective semi- 
classical phase-space flow obtained by Peierls substitution, and involve a slowly 
varying envelope whose dynamics is governed by a homogenized Schrodinger 
equation with time-dependent effective mass. The corresponding adiabatic de- 
coupling of the slow and fast degrees of freedom is shown to be valid up to 
Ehrcnfcst time scales. 



1. Introduction 

1.1. General setting. We consider the following semiclassically scaled Schrodinger 
equation: 

{2 
ied t r + £ -M E = V T (I) r + V{xW, (t,x) elx R d , 
^jt=o = $>> 

with d ^ 1, the spatial dimension, and ip E = ip £ (t,x) € C. Here, we already have 
rescaled all physical parameters such that only one semiclassical parameter e > 
(i.e. the scaled Planck's constant) remains. In the following we shall be interested 
in the asymptotic description of ip £ (t,x) for e -C 1. To this end, the potential 
Vr(y) 6 t is assumed to be smooth and periodic with respect to some regular 
lattice T ~ Z d , generated by a given basis {771, . . . , rja}, rjt € R d , i.e. 

(1-2) V T (y + 1 ) = V r (y), VyeR d , 7 e.r 

where 



a 

{7 = ^7^eM d : 7 fez}. 



r = i - = 

In addition, the slowly- varying potential V is assumed to satisfy the following: 

Assumption 1.1. V G C 3 (M d ;M) and d^V e L°° (R d ) , for | 7 | = 2,3. 

Note that this implies that V(x) grows at most quadratically at infinity. Equa- 
tion [TTT] describes the dynamics of quantum particles in a periodic lattice-potential 
Vr under the influence of an external, slowly varying driving force F = — VV(x). 
A typical application arises in solid state physics where (|1.1[) describes the time- 
evolution of electrons moving in a crystalline lattice (generated by the ionic cores) . 
The asymptotics of (jl.ip as e — > 0+ is a natural two-scale problem which is well- 
studied in the physics and mathematics literature. Early mathematical results are 
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based on time-dependent WKB type expansions [3J Q]3 [37J (see also [7J for a more 
recent application in the nonlinear case), which, however, suffer from the appear- 
ance of caustics and are thus only valid for small times. In order to overcome this 
problem, other methods based on, e.g., Gaussian beams [11], or Wigner measures 
[131 [TJ], have been developed. These approaches yield an asymptotic description 
for time-scales of order 0(1) (i.e. beyond caustics). More recently, so-called space- 
adiabatic perturbation theory has been used (together with Weyl pseudo-differential 
calculus) to derive an effective Hamiltonian, governing the dynamics of particles in 
periodic potentials Vr under the additional influence of slowly varying perturbations 
[221 135] . The semi-classical asymptotics of this effective model is then obtained in 
a second step, invoking an Egorov-type theorem. 

On the other hand, it is well known that in the case without periodic potential, 
semiclassical approximations which are valid up to Ehrenfest time t ~ 0(lnl/£) 
can be constructed in a rather simple way. The corresponding asymptotic method 
is based on propagating semiclassical wave packets, or coherent states, i.e. ap- 
proximate solutions of (|1.1[) which are sufficiently concentrated in space and in 
frequency around the classical Hamiltonian phase-space flow. More precisely, one 
considers 

(1.3) ^(i, x) « e- d ^u (t, ^-§^) e <(S(t)+(x-,(t)).p(t))/ e) 

where (q(t),p(t)) € C 3 (R;R 2d ) satisfies Hamilton's equation of motion 



(1.4) 



q(t)=p(t), q(0) = q O , 

p(t) = -V x V(q(t)), p(0)=po, 



and the purely time-dependent function S(t) denotes the classical action (see H1.3\ 
below). The right hand side of (jl.3p corresponds to a wave function which is equally 
localized in space and in frequency (at scale \fs), so the uncertainty principle is 
optimized. In other words, the three quantities 

-?(*). 



, and 

L 2 (R d ) 



have the same order of magnitude, 0(1), as e — > 0. The basic idea for this type 
of asymptotic method can be found in the classical works of [16l [27] (see also 
[H [29] for a broader introduction). It has been developed further in, e.g., [5J [HI 
[311 1351 135] and in addition also proved to be applicable in the case of nonlinear 
Schrodinger equations [B] (a situation in which the use of Wigner measures of space- 
adiabatic perturbation theory fails) . Asymptotic results based on such semiclassical 
wave packets also have the advantage of giving a rather clear connection between 
quantum mechanics and classical particle dynamics and are thus frequently used in 
numerical simulations (see e.g. [T2"]b 

Ehrenfest time is the largest time up to which the wave packet approximation is 
valid, in general. Without any extra geometric assumption, the coherent structure 
may be lost at some time of order C In 1/e, if C is too large. See e.g. [5J [TUl [30J [35J 
135] and references therein. 

Interestingly enough, though, it seems that so far this method has not been 
extended to include also highly oscillatory periodic potentials Vr (f), and it will be 
the main task of this work to do so. To this end, it will be necessary to understand 
the influence of Vr (§) on the dispersive properties of the solution tp £ (t,x). In 
particular, having in mind the results quoted above, one expects that in this case 
the usual kinetic energy of a particle E = \\k\ 2 has to be replaced by E m (k), i.e. 
the energy of the m-th Bloch band associated to Vr- In physics this is known under 
the name Peierls substitution. We shall show that under the additional influence 
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of a slowly varying potential V(x), this procedure is in fact asymptotically correct 
(i.e. for e < 1) up to Ehrenfest time, provided the initial data tpQ is sufficiently 
concentrated around (qa,Po) G K 2d - 

Remark 1.2. Indeed, we could also allow for time-dependent external potentials 
V(t, x) G K measurable in time, smooth in x, and satisfying 

d2V £L°° (R t xR d x ) : |7|=2,3. 

Under this assumptions, it is straightforward to adapt the analysis given below. 
For the sake of notation, we shall not do so here, but rather leave the details to the 
reader. 

1.2. Bloch and semiclassical wave packets. In order to state our result more 
precisely, we first recall some well-known results on the spectral theory for periodic 
Schrodinger operators, cf. [5T1 HO] : 

#per :=-~A y + V r {y). 

Denote by Y C T the (centered) fundamental domain of the lattice T, equipped 
with periodic boundary conditions, i.e. Y ~ T d . Similarly, we denote by Y* ~ T d 
the fundamental domain of the corresponding dual lattice. The latter is usually 
referred to as the Brillouin zone. Bloch-Floquet theory asserts that H per admits a 
fiber-decomposition 

H p™ = J H r{k) dk, 
where for k G Y*, we denote 

H r (k) = i (-*V„ + k) 2 + V r (y) . 
It therefore suffices to consider the following spectral problem on Y: 
(1.5) H r (k) Xm (;k) = E m (k) Xm {; k), k e Y*, m e N, 

where E m (k) 6 R and Xm{y,k), respectively, denote an eigenvalue/eigenvector 
pair of Hr(k), parametrized by k € Y* , the so-called crystal momentum. These 
eigenvalues can be ordered increasingly, such that 

Ex{k) < .. . < E m {k) < E m+1 (k) < . . . , 

where each eigenvalue is repeated according to its multiplicity (which is known to 
be finite). This implies that 

specif per = |J {E m (k) ; k e Y*} C K, 

where {E m (k); e F} is called the m-th energy band (or Bloch band). The 
associated eigenfunctions Xmiy^k) are r*-periodic w.r.t. k and form a complete 
orthonormal basis of L 2 (Y). Moreover, the functions Xm('>k) € H 2 (Y) are known 
to be real-analytic with respect to k on y*\fi, where is a set of Lebesgue measure 
zero (the set of band crossings). 

Next, we consider for some to S N the corresponding semi- classical band Hamil- 
tonian, obtained by Peierl's substitution, i.e. 

h%(k,x) = E m (k) + V(x), (k,x) g Y* x R d , 

and denote the semiclassical phase space trajectories associated to by 

f - V fe S ro (p(i)) , g(0) = go , 

( ' } \*(*)=-V»V (?(*)), p(0)=po. 

This system is the analog of (|1.4p in the presence of an additional periodic potential. 
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Example 1.3 (No external potential). In the case V(x) — 0, we simply have 

(1.7) P(t)=p , q(t) = q +tV k E m (po), 
that is, a shift with constant speed lo = V E m (po). 

In order to make sure that the system (|1.6[) is well-defined, we shall from now 
on impose the following condition on E m (k). 

Assumption 1.4. We assume that E m (p(t)) is a simple eigenvalue, uniformly for 
all t£R, i.e. there exists a 6 > 0, such that 

\E rn (p(t)) - E n (k)\ > *, V^m,teR,fc€ Y*. 

It is known that if E m (k) is simple, it is infinitely diffcrentiable and thus the right 
hand side of (ll.6[) is well defined. Under Assumption II .41 we consequently obtain a 
smooth semi-classical flow (qo,po) i-> (q(t),p(t)), for all t € R. In addition, one can 
choose Xm{y, k) to be T-periodic with respect to y and such that (y, t) i-> XmiViPif)) 
is bounded together with all its derivatives. 

Example 1.5. By compactness of Y* , Assumption II .41 is satisfied in either of the 
following two cases: 

(i) If E m (k) is a simple eigenvalue for all k £ Y* . In particular, in d = 1 it is 
known that every E m (fc) is simple, except possibly at k — or at the edge 
of the Brillouin zone. 

(ii) If V(x) = and E m (k) is simple in a neighborhood of k = po (which is 
sufficient in view of Example 1 1 . 3[) . 

1.3. Main result. With the above definitions at hand, we are now able to state 
our main mathematical result. To this end, we first define a semiclassical wave 
packet in the m-th Bloch band (satisfying Assumption II .4[) by 

(1.8) <f(t,x) = e-^u (t, ^§^) Xm (*p(t)) e*»<*->/« 

with q(t),p{t) given by system (jl.6p and u{t,z) £ C, a smooth slowly varying 
envelope which will determined by an envelope equation yet to be derived (see 
below). In addition, the e-oscillatory phase is 

(1.9) m (i, x) = S m (t) + p(t) ■ (x ~ q(t)), 

where S m (t) £ R is the (purely time-dependent) semi-classical action 

(1.10) S m (t) = f L m (p( S ),q( S ))ds, 

Jo 

with L rn denoting the Lagrangian associated to the effective Hamiltonian h^, i.e. 

L m (p(s), q(s)) = p{s) ■ q(s) - h% (p(a), q(s)) 
( ' } =p(s)-VE m (p(s))-h^(p(s),q(s)), 
in view of (|1.6[) . 

Remark 1.6. Note that this is nothing but the Legendre transform of the effective 
Hamiltonian h^. As in classical mechanics, one associates to a given Hamiltonian 
H(p, q) a Lagrangian via L(p, q) = p ■ q — H(p, q). 

The function ip £ given by (ll.8[) generalizes the usual class of semiclassical wave 
packets considered in e.g. [THl [37]. Note that in contrast to two-scale WKB ap- 
proximation considered in [3j [15j [37] , it involves an additional scale of the order 
l/\/s, the scale of concentration of the amplitude u. In addition, (|1.8|) does not 
suffer from the appearance of caustics. Nevertheless, in comparison to the highly 
oscillatory Bloch function Xm , the amplitude is still slowly varying and thus we can 
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expect an adiabatic decoupling between the slow and fast scales to hold on (long) 
macroscopic time-scales. Indeed, we shall prove the following result: 

Theorem 1.7. Let Vr be smooth and V satisfy Assumptions ["PI In addition, let 
Assumption \l .4\ hold and the initial data be given by 

with qo,po £ M. d and some given profile uo £ S(M. d ). Then there exists C > such 
that the solution of (jl.ip can be approximated by 

y/(t)~^(t)\\ L2(Rd) ^C^e ct . 

Here, (p e is given by (|1.8I) with 

u(t, z) — v(t, z) exp f3(r)dT^j , 

where f3(t) £ iR is the so-called Berry phase term 

P(t) := (Xm(p(t)),V kX m(j>(t)))mY) ■ W(q(t)), 
and v £ C(R; <5>(R d )) satisfies the following homogenized Schrodinger equation 

(1.12) id t v+~dW z (VlE m {p{t))-V z )v=^(z,VlV(q(t))z)v, v\ t=Q = u . 
In particular there exists Cq > so that 

sup ||^ e (*)-<p e (*)||i,3(R<*) — >0. 

O^t^Co In J 

Remark 1.8. In fact it is possible to prove the same result under less restrictive 
regularity assumptions on uq and Vr . Indeed, Proposition 15.11 shows that it is suf- 
ficient to require that uq belongs to a certain weighted Sobolev space. Concerning 
the periodic potential, it is possible to lower the regularity considerably, depending 
on the dimension. For example, in d = 3 it is sufficient to assume Vr to be infinitcs- 
imally bounded with respect to -A. This implies x m (-,fc) £ H 2 (R 3 ) <-> L°°(]R 3 ), 
which, together with several density arguments (to be invoked at different stages 
of the formal expansion), is enough to justify the analysis given below. 

Theorem 11.71 provides an approximate description of the solution to (jl.ip up to 
Ehrenfest time and can be seen as the analog of the results given in [TBI HZ1 HI 
IH1 133 133 [3E] where the case of slowly varying potentials V(x) is considered. 
The proof does not rely on the use of pseudo-differential calculus or space space- 
adiabatic perturbation theory and can thus be considered to be considerably simpler 
from a mathematical point of view. In fact, our approach is similar to the one 
given in [TB], which derives an analogous result for the so-called Born-Oppenheimer 
approximation of molecular dynamics. Note however, that we allow for more general 
initial amplitudes, not necessarily Gaussian. Indeed, in the special case where the 
initial envelope uq is a Gaussian, then its evolution u remains Gaussian, and can 
be completely characterized; see ^4.31 Also note that in contrast to the closely 
related method of Gaussian beams presented in, e.g., [TT], we do not need to include 
complex-valued phases and in addition, obtain an approximation valid for longer 
times. 

The Berry phase term is an example for so-called geometric phases in quantum 
mechanics. It is a well known feature of semiclassical approximation in periodic 
potentials, see, e.g., [3S] for more details and a geometric interpretation. The 
homogenized Schrodinger equation features a rather unusual dispersive behavior 
described by a time-dependent effective mass tensor M(t) — V\E m (p(t)), i.e. the 
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Hessian of E m (k) evaluated at k = p(i). To our knowledge, Theorem II .71 is the first 
result in which a Schrodinger type equation with time- dependent effective mass has 
been rigorously derived (see also the discussion in Remark 13.11) . 

Remark 1.9. Let us also mention that the same class of initial data has been con- 
sidered in [T] for a Schrodinger equation with locally periodic potential Vr(x, y) and 
corresponding ^-dependent Bloch bands E m (k; x). In this work, the authors derive 
a homogenized Schrodinger equation, provided that ipQ is concentrated around a 
stationary point point xq,pq of the semiclassical phase space flow, i.e. 

VkE m (p ;q ) = V x E m (p ;q ) = 0. 

This implies q(t) — qo and p(t) — po, for all t g R, yielding (at least asymptotically) 
a localization of the wave function. We observe the same phenomenon in our case 
under the condition V(x) — and VkE m (k) = (see Example 1 1 .3[) . 

This work is now organized as follows: In the next section, we shall formally 
derive an approximate solution to (11.11) by means of a (formal) multi-scale expan- 
sion. This expansion yields a system of three linear equations, which we shall solve 
in Section [3] In particular, we shall obtain from it the homogenized Schrodinger 
equation. The corresponding Cauchy problem is then analyzed in Section 21 where 
we also include a brief discussion on the particularly important case of Gaussian 
profiles (yielding a direct connection to 16 ). A rigorous stability result for our 
approximation, up to Ehrenfest time, is then given in Section [5j 

Remark 1.10. We expect that our results can be generalized to the case of 
(weakly) nonlinear Schrodinger equations (as considered in [SJ [7]). This will be 
the aim of a future work. 



2. Formal derivation of an approximate solution 

2.1. Reduction through exact computations. We seek the solution ij) e of (jl.ip 

in the following form 

(2.1) ^{t,x)=e- d ' i W ^,^-|^,^ e ^('^, 

where the phase 4>(t,x) is given by (|1.9I) . the function U e ~ U e (t,z,y) is assumed 
to be smooth, T-periodic with respect to y, and admits an asymptotic expansion 

(2.2) U%t,z,y) ~ yy^i.z.j/). 

Note that due to the inclusion of the factor e _d / 4 the L 2 (M. d ) norm of the right hand 
side of (12.11) is in fact uniformly bounded with respect to e, whereas the L°°(M. d ) 
norm in general will grow as e — > 0. The asymptotic expansion 12 . 2 1 therefore has to 
be understood in the L 2 sense. 

Taking into account that in view of (jl.9l) . V ' x 4> m (t,x) = p(t), we compute: 

£ d/* e -i<t> m / He g t Tpe = i£ g tU e _ f ^ . ^ jje _ d t (j) m U E , 
£ d/4 e -*4> m /s £ 2 A ^e = £AzW + AyU e + 2 ^ {V y . V ' Z )W - \p\ 2 W 

+ 2i^/ep ■ V z U e + Zip ■ V y W, 

where in all of the above expressions, the various functions have to be understood 
to be evaluated as follows: 

t// = ^(t,x) ; W=W (t *~ q ^\ - 
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2- 



b E = -dt<kmU e + -A y W - -\p\ 2 U £ + ip ■ V V U £ - V T (y)U £ - V(q)W 



Thus, ordering equal powers of e in equation (|1.9p we find that 

£ d/i e - l(t > m /e L edt ^ + _ Vr Q ^ - V(x)^ = b £ a + y/etil + eb 

with 

b\ = -iq ■ V Z U £ + (V„ • VzW + ip ■ V 2 W e , 

b £ 2 = id t W + \^zU E . 

So far, we have neither used the fact that q(t), p(t) are given by the Hamiltonian 
flow p. 61) . nor the explicit dependence of <j> m on time. Using these properties, allows 
us to rewrite 6g, b\ , b\ as follows: 

b% = (K c n ( P (t),q(t)) + W (q(t)) ■ (x - q{t)))W - H T (p{t))W - V(q(t))U s , 
b\ = i (p(t) - V k E m (p(t))) • V z W e + (V„ • V z )W e , 

Now, recall that in the above lines, U e is evaluated at the shifted spatial variable 
z = (x — q(t))/*/s. Taking this into account, we notice that the above hierarchy has 
to be modified, and we find: 

b £ o = h%(p(t),q(t))U e - H r { P {t)W - V (q(t) + zji) U £ , 

b\ = i (p(t) - V k E m (p(t))) • V Z U £ + (V„ ■ V Z )U £ + (W (q(t)) ■ z)U\ 

b\ = id t w + \&zU £ . 

Next, we perform a Taylor expansion of V around the point q(t): 

V (q(t) + z^)=V (q(t)) + ViW (q(t)) -z+ £ -(z, VV (q(t)) z) + O (e 3 / 2 (zf) , 

since V is at most quadratic in view of Assumption ll.il Recalling that h^(p, q) — 
E m (p) + V(q), the terms involving V(q) cancel out in 6q, the terms involving W(q) 
cancel out in b\, and thus, we finally obtain: 



Lemma 2.1. Let the Assumptions [777J \l-4\ hold and i\j £ be related to IA £ through 
AO). Then it holds 



2 

ied t i) £ + £ —^ £ - V v ( j) ijf - V(x)iJj £ = 

(bl + ^ebl+ebl + e^ 2 ^) (t,z,y) 



e i<j>m./e 

~7m 



(z,y)=(- 



with 



b £ ^(E m (p(t))-H r (p(t)))U £ 1 

b\ = i (p(t) - V k E m (p(t))) • V Z U £ + {V y ■ V Z )U £ , 

b £ 2 = id t U £ + ^A Z U £ - i (z, V 2 V (q(t)) z) U £ , 
and a remainder r £ (t,z,y) satisfying 

\r £ (t,z,y)\^C(z) 3 \U £ (t,z,y)\, V(t, z, y) e R x R d x Y, 
where the constant C > is independent oft, z, y and e. 
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2.2. Introducing the approximate solution. We now expand W in powers of 
e, according to (|2.2p . To this end, we introduce the following (time-dependent) 
linear operators 

L = E m (p(t)) - H r (p(t)) , 

Li = i (p(t) - V k E m (p(t))) • V z + V„ • V z , 

i 2 =*a t + iA z -i(^V 2 V(<z(t))z). 

In order to solve (jl.ip up to a sufficiently small error term (in L 2 ), we need to 
cancel the first three terms in our asymptotic expansion. This yields, the following 
system of equations 

f L Q Uo = 0, 
(2.3) J 

[ L Q U 2 + + L 2 U Q = 0. 

Assuming for the moment that we can do so, this means that we (formally) solve 
(II. ip up to errors of order e 3 / 2 (in L 2 ), which is expected to generate a small 
perturbation of the exact solution (in view of the e in front of the time derivative 
of V> £ in (TTTTj) ). 

We consequently define the approximate solution 



(2 



.4) ^ pp (f, x) := e^ 4 (C/o + y/iUi + eU 2 ) (t, ^ e *~(t,«)/« 



In view of Lemma l2.1[ we thus have the following result (provided we can solve the 
system (|2.3[) in a unique way): 

Lemma 2.2. Let ifil pp given by (|2.4j) . where Uq,U\,U 2 solve (|2.3[) . T/ien 

(fe$ + y A - y r - ^ pp = e3/2 ^ + r t + r I) (*, *, y) 

where the remainder terms rf , r| are given by 

rl(t,z,y) = L 2 Ui(t,z,y), r%{t,z,y) = LiU 2 (t,z,y), 
and r E satisfies 

\r £ (t, z,y)\^C (z) 3 | (C/ + yfeU x + eU 2 ) (t, z, y) | , V(t, z, y) G R x M d x F, 
where the constant C > is independent of t,z,y and e. 

3. Derivation of the homogenized equation 

3.1. Some useful algebraic identities. Given the form of Lo, the equation 
LqUq — implies 

(3.1) U (t,z,y) =u(t,z)xm(y,p(t)) ■ 

Before studying the other two equations, we shall recall some algebraic formulas 
related to the eigenvalues and eigenvectors of Hp. First, in view of the identity 
(|1.5[) . we have 

(3.2) V fc (H r - E m ) Xm + (H r - E m ) V fcX ™ = 0. 
Taking the scalar product in L 2 (Y) with Xm, we infer 

Vfc-E m = (Xmi VfciJrXm) L 2 (Y) + (Xmi (#r ~ -Em) VfeX m ) £2( Y ) • 

Since iJr is self-adjoint, the last term is zero, thanks to (jl.5D . Wc infer 

(3.3) V k E m (k) = (xm, (-iVj, + fc) Xm) £ 2 (y) ■ 



SEMICLASSICAL WAVE PACKETS IN PERIODIC POTENTIALS 9 

Differentiating (I3.2j) again, we have, for all j, £ E {1, . . . , d}: 

dljke — E m ) Xm + dkj {Hr — E m ) dk e Xm + dk t {Hr — E m ) d^-x-m 
+ (H r -E m )d 2 k]keXm = 0. 
Taking the scalar product with % m , we have: 

d Lk e E m(k) = Sji + {{-id Vj + kj) d kl; Xm + {-id ye + h) d k] Xm, Xm) L 2 (Y) 

( 3 - 4 ) 

- (dkeEmdkjXm + dk 3 E m dk e Xm, Xm) l2( y) ' 

3.2. Higher order solvability conditions. By Fredholm's alternative, a neces- 
sary and sufficient condition to solve the equation LqUi + L\Uq — 0, is that L\Uq 
is orthogonal to kerLo) that is: 

(3-5) (Xm,LiU ) L 2(Y) = 0. 

Given the expression of L\ and the formula (|3.1[) . we compute 

L x u a = % (p(t) - V k E m (p(f))) • v«tt(t, z) Xm (y,p(t)) + v yXm {y,p{t)) ■ V z u{t, z). 

In view of (|3.3[) . we infer that (|3.5[) is automatically fulfilled. We thus obtain 

Ui(t,z,y) = Ui(t,z)xm (y,p{t)) + Ui(t,z,y), 

where u^, the part of XJ\ which is orthogonal to kerLo> is obtained by inverting an 
elliptic equation: 

ui = -L^LxUq. 
Note that the formula for LiUq can also be written as 

LxC/q = -*V fe {E m (p(t)) - H r (p(t))) X m (y,p(t)) • V z u(i, z), 
thus taking into account (|3.2[) . we simply have: 

ui(t,z,y) = -iV k Xm (y,p(t)) ■ V z u(t,z). 

At this stage, we shall, for simplicity choose u\ = 0, in which case U\ becomes 
simply a function of u: 

(3.6) U 1 (t,z,y) = -W k \ in {y,p{t))-V z u{t,z). 
As a next step in the formal analysis, we must solve 

L U 2 + LiUi + L 2 U = 0. 
By the same argument as before, we require 

(3.7) {Xm,L 1 U 1 +L 2 U () ) L2(Y) ={). 
With the expression (|3.6[) . we compute 

d 

LiUi = {(Pi*) - VkE m (p(t))) j d ke xm (y,p(t)) - idi. kt Xm (y,p(t))) d 2 Z]Ze u, 

and we also have 

L 2 U = ( (id t + \^--\ <*> V 2 F (q(t)) z)) u) Xm (y,p(t)) 
+ iup(t) • Vfex m (y,p(t)). 
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Recalling that p(t) — — W (?(£)), we find: 

(Xm,LiUi + L 2 Uq) l 2 {Y) = 

= (id t + i A, - i (z, V V z)^ u - *W(?(t)) • (x™, V feXm ) L2(y) u 

- X! (Xm, <9 fej £ m (p(t)) <9fc,Xm + id k 3 k e X m ) d 2 Z]Ze U 

By making the last sum symmetric with respect to j and £, and using (|3.4I) . we 
finally obtain the homogenized Schrodinger equation with time-dependent effective 
mass: 

(3.8) id t u + idiv, (V 2 k E m (p(t)) • V.) u = ^ («, V 2 V ( ? (t)) z) u + if3(t)u, 
where we recall that 

(3(t) = W(g(i)) ■ (x TO (p(i)), V fe Xm(p(i))) i2( r) - 
the so-called Berry phase term. From ||Xm||L 2 (Y) = 1, we infer that 

Re(Xm, VfcXm) L 2 (y) = 0. 

In other words, (Xm,VtXm)i2(y) G iM and thus i/3(t) G K, acts like a purely time- 
dependent, real-valued, potential. Thus, invoking the unitary change of variable 

v(t,z) = u(t, z)e~ P^ ds 

implies that v(t,z) solves (|1.12[) . Equation (13.81) models a quantum mechanical 
time-dependent harmonic oscillator, in which the time dependence is present both 
in the differential operator, and in the potential. 

Remark 3.1. In the case where V(x) — 0, we have /3(t) = and pit) = po (in 
view of Example II .3[) . In this case v(t, z) = u(t, z) and Equation (|3.8[) simplifies to 
an equation with a time-independent effective mass tensor: 

id t u + idiv z (V 2 k E m (p ) ■ V,) u = 0. 

This equation has been derived in [5] (see also [3111 [23 US] for similar results). 
Note, however, that in the quoted works the scaling of the original equation (|l.ip 
is different (i.e. not in semiclassical form). 

Assuming for the moment that we can solve (|3.8p . we can write 

(3.9) U 2 {t,z,y) = u 2 {t,z)x m (y,p(t)) +u 2 L (t,z,y), 
where 

«a = - L o' i L iUi + L 2Uo) ■ 

Like we did for we shall from now on also impose u 2 = and thus U 2 = u%. 
For the upcoming analysis, we define the following class of energy spaces 

S fe - | / G L 2 (R d ) ; := £ ||^^/|| L2(Rd) < oo; fc e N 

Having in mind f3TT]l , (|3~o) . (j3~^|) . and the fact that L 1 : L 2 (Y) -> iJ 2 (F), we 
directly obtain the following result: 
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Proposition 3.2. Let Assumption \l .4\ hold and let u £ C(R;£ fe ) solve (|3.8[) . Set 

Uo(t, z, y) = u{t, z)x m (y,p{t)) , 
Ui(t,z,y) = -iVkXm, (y,p(t)) ■ V z u(t,z), 
U 2 (t, z, y) = -L 1 (Lil7i(t, z, y) + L 2 U (t, z, y)) . 
Then Uj £ C(R; x W°°'°°{Y)), for j = 0, 1, 2 and (Z7 0) C7i, C/ 2 ) solves (1231) . 

4. The envelope equation 

We examine the Cauchy problem for Q3.8p . with special emphasis on the large 
time control of u. 

4.1. The general Cauchy problem. Equation 13.81 can be seen as the quantum 
mechanical evolutionary problem corresponding to the following time-dependent 
Hamiltonian, 

(4.1) H(t, z, = ~ (C, V\E m (pit)) C) + \ (zV 2 V (q(t)) z) + i8(t). 

Under Assumptions ITTTI and 1X741 this Hamiltonian is self-adjoint, smooth in time, 
and quadratic in (z,Q (in fact, at most quadratic would be sufficient). Using the 
result given in [5H p. 197] (see also [IS]), we directly infer the following existence 
result: 

Lemma 4.1 (From [24]). For d 1 and vo G L 2 (R d ), consider the equation 

(4.2) iftv+i a jk(t)djk v= ^ b jk {t)xjX k \ ; V| t=0 = v . 

l<j,fe<d l^j.fcsjd 

If the coefficients ajk and bjk are continuous and real-valued, such that the matrices 
( a jk)j,k and (bjk)j,k are symmetric for all time, then (14.21) has a unique solution 
v € C(R; L 2 (R d )).' It satisfies 

(4.3) \\v(t)\\ L 2 m = llvoll^cRd), VteM. 

Moreover, if vq £ Y, k for some k £ N, then v £ C(R;I] fc ). 

In particular, this implies that if uo £ then (|1.12p has a unique solution 
v £ C(R;S fc ). As a consequence, (|3.8[) has a unique solution w G C(IR;E fe ) such 
that u\t—Q = uq. 

Remark 4.2. It may happen that the functions ajk are zero on some non-negligible 
set. In this case, (|4.2p ceases to be dispersive. Note that the standard harmonic 
oscillator is dispersive, locally in time only, since it has eigenvalues. We shall see 
that this is not a problem in our analysis though. 

4.2. Exponential control of the envelope equation. To prove Theorem 11.71 
we need to control the error present in Lemma 12.21 for large time. In general, i.e. 
without extra geometric assumptions on the wave packet, exponential growth in 
time must be expected: 

Proposition 4.3. Let £ E fe , k £ N. Then the solution u to (|3.8p satisfies 
u £ C(R; S fc ), and there exists C > such that 

\\u(t,-)\\ sh < Ce c \ t^O. 

Proof. The result can be established by induction on k. The constant C must 
actually be expected to depend on k, as shown by the case of 

1. \z\ 2 
io t u+-A z u = — — u ; U|f = o = uo. 
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There, the fundamental solution is explicit (generalized Mehler formula, see e.g. 
[21]), and we check that ||«(t)|| E j« behaves like e kt . 

For k = 0, the result is obvious, since in view of Lemma 14.11 the L 2 -norm is 
conserved. The case k = 1 illustrates the general mechanism of the proof, and 
we shall stick to this case for simplicity. The key remark is that even though 
the operators z and V z (involved in the definition of S 1 ) do not commute with 
the Hamiltonian (I4.1[) . the commutators yield a closed system of estimates. First, 
multiplying (I3.8[) by z, we find 

(id t -H)zu = - [H, z]u = V 2 E m (p(i)) \7 z u. 

The L 2 estimate (|4.3I) then yields 

||;zu(t)|| i2 < H^uollia + / \\\7 2 k E m (p(s))\7 z u(s)\\ L2 ds 



< ||zu || L 2 + C j \\V z u(s)\\ L2 ds, 
Jo 



for some C independent of t, since ^J\E m is bounded on Y* which is compact. 
Similarly, 

(id t ~ H) V z u = - [H, V z ] u = V 2 V (q(t)) zu, 
and, in view of Assumption 11.11 



\\V z u(t)\\ L2 < ||V z it || L2 + C / ||jMi(«)|| za ds. 

Jo 

Summing over the two inequalities and using the conservation of mass, we infer 
II«(*)IIei < II «Q II E i + c \\ u ( s )hi ds, 



and Gronwall's lemma yields the proposition in the case k — 1. By induction, 
applying (z, V z ) to p.8p k times, the defects of commutation always yield the same 
sort of estimate, and the proposition follows easily. □ 



4.3. Gaussian wave packets. In the case where the initial datum in (|3 .8[) is a 

Gaussian, we can compute its evolution and show that it remains Gaussian, by 
following the same strategy as in [TB] (see also J71 US]). As a matter of fact, the 
order in which we have proceeded is different from the one in |16) . since we have 
isolated the envelope equation (|3.8[) before considering special initial data. As a 
consequence, we have fewer unknowns. Consider (|3.8I) with initial datum 



(4.4) "(°^)=(d^ eX p(4^^))' 
where the matrices A and B satisfy the following properties: 

(4.5) A and B are invertiblc; 

(4.6) BA~ X is symmetric : BA~ X = Mi + 1M2, with Mj real symmetric; 

(4.7) KeBA^ 1 is strictly positive definite; 

(4.8) (ReBA^ 1 ) 1 = AA* . 

Proposition 4.4. Let u solve (|3.8p . with initial datum (|4.4p . where the matrices 
A and B satisfy (|4.5[) - (|4.8[) . Then for all time, u(t,z) is given by 

(4.9) „(*,,) - (det ^ )) i /2 exp {-\ (z,B(t)A(t)-^ , 
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where the matrices A(t) and B(t) evolve according to the differential equations 
{A{t)=iV 2 E rn {p{t))B(t) ■ A(0)=A, 
\B(t)=iVlV(q(t))A(t) ; B(0) = B. 

In addition, for all time teR, A(t) and B(t) satisfy (|4.5[) - (|4.8j) . 

Proof. The argument is the same as in [IB] (see also jTTJ \TE\ ) : One easily checks that 
if A(t) and B(t) evolve according to (|4.10p . then u given by (|4.9[) solves (I3.8[) . On the 
other hand, it is clear that (|4.10l) has a global solution. Finally, since Vf.E m (p(t)) 
and ^7 2 V (q(t)) are symmetric matrices, it follows from [THl Lemma 2.1] that for all 
time, A(t) and B{t) satisfy (g3]) -g^]). □ 

5. Stability of the approximation up to Ehrenfest time 

As a final step we need to show that the derived approximation V'appW indeed 
approximates the exact solution ip E (t) U P to Ehrenfest time. 

Proposition 5.1. Let Assumptions QT7] and \l-4\ hold and uq € S 5 . Then there 
exists C > such that 

n E (t)-^ pp (t)\\ L2(Rd) ^CV^e ct , 

where V'appfA x) is given by (|2.4[) with u(t,z) solving (|3.8p subject to u\ t=Q = uq, 
and Uq, U\, XJi are given by Proposition[ 



Proof. First, note that ip e and ?Pl pp do not coincide at time t = 0, since elliptic 
inversion has forced us to introduce U\ and U2 which are not zero initially. Setting 
w e — 4> £ — ipl pp , and using and Lemma \2. 21 we see that the error solves 



[ied t + yA — Vr — V S jw e = ~^> 2 (r E +r{+ rf) (t, z, „) (z y)=( ^, jf } . 

™ £ (0, x) = e' d / 4 (y/IUi + eU 2 ) (0, ^B, |) (».*>/«. 

The a-priori L 2 estimate yields 

IKWIIl^ < V^\\Ui(0)h*L~ +e\\U 2 (0)\\ LlLT 

+ V~ef (||r £ ( S )|| w + \\r\{s)\\ LlLT + K(a)|| zj£ .) 



ds. 

The assertion then follows from Lemma 12.21 (establishing the needed properties 
for the functions r e , rf and r^), Proposition 13.21 and Proposition 14.31 With this 
approach, we need to know that r £ is in L 2 , so Uq, U\ and U2 have three momenta 
in L 2 : in view of Proposition 13.21 and Proposition 14.31 this amounts to demanding 

u e s 5 . □ 

This asymptotic stability result directly yields the assertion of Theorem 11.71 
End of the proof of Theorem \1.7\ To conclude, it suffices to notice that 

ipFQt, x) = yj £ app (t, x) - (ViC/i + eU 2 ) (t, z, y) 

so Proposition 13.21 and Proposition [43] imply 



(*,y)=( £ 



k e (*)-Vapp(i)IU=><<Vie 



ct 



This estimate and Proposition 15. II yield Theorem 11.71 □ 
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Remark 5.2. The construction of the approximate solution ^>| pp has forced us 
to introduce non-zero correctors U± and U2, given by elliptic inversion. Therefore, 
we had to consider well-prepared initial data for ipl pp . This aspect is harmless 
as long as one is interested only in the leading order behavior of ip £ as e — > 0. 
As a consequence, our approach would not allow us to construct arbitrary accurate 
approximations for ip e (in terms of powers of e), unless well-prepared initial data are 
considered, i.e. data lying in so-called super- adiabatic subspaces, in the terminology 
of [25] (after [H]). This is due to the spectral analysis implied by the presence of 
the periodic potential Vr, and shows a sharp contrast with the case Vr = 0. 

Of course the above given stability result immediately generalizes to situations 
where, instead of a single ip E , a superposition of finitely many semiclassical wave 
packets is considered, 

n—l V V / 

Since the underlying semiclassical Schrddinger equation is linear, each of these 
initial wave packets will evolve individually from the rest, as in Theorem 11.71 Up 
to some technical modifications, it should be possible to consider even a continuous 
superposition of wave packets, yielding a semiclassical approximation known under 
the name "frozen Gaussians" , see [19] . 
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